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We investigate the ground state properties of an equal mixture of two species of bosons in its 
Mott-insulator phase at a filUng factor two per site. We identify one type of spin triplet-singlet 
transition through the competition of ground state. When the on-site interaction is weak {U < Uc) 
the two particles prefer to stay in the lowest band and with weak tunnelling between neighboring 
sites the system is mapped into an effective spin-1 ferromagnetic exchange Hamiltonian. When the 
interaction is tuned by a Feshbach resonance to be large enough (U > (7c), higher band will be 
populated. Due to the orbital coupling term S'^S~ in the Hamiltonian, the two atoms in different 
orbits on a site would form an on-site singlet. For a non-S'[/(2)-symmetric model, easy-axis or 
easy-plane ferromagnetic spin exchange models may be realized corresponding to phase separation 
or counter-flow superfluidity, respectively. 

PACS numbers: 05.30.Jp, 03.75.Mn, 73.43.Nq 
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The study of quantum phase transition in optical lat- 
tices has made great progress both theoretically and ex- 
perimentally and becomes one of focusing issues of 
current interest in the exploration of rich physics in ul- 
tracold atomic systems. Jaksch et. al. predicted that the 
dynamics of a single-component Bose gas loaded into the 
lowest band of an optical lattice is well described by the 
Bose-Hubbard model [l| and Greiner et. al. experimen- 
tally confirmed that the phase transition from superfluid 
phase to Mott-insulator could be realized by suppress- 
ing tunnelling between neighboring sites [2]. For single- 
component bosons without internal degrees of freedom, 
the superfluid-insulator transition in a periodic lattice 
has been extensively studied by various methods [H, H, 01 ■ 
When the spinless bosons are in the Mott phase, the 
on-site fluctuation of particle numbers is suppressed [3]. 
Many studies have shown that multi-component bosonic 
or fermionic gases in optical lattices exhibit much richer 
phase diagrams 

BIS Mill- An intri guing feature 
of the multi-component Bose systems is the structure of 
their internal "spin" degree of freedom. The recent for- 
mation of bound repulsive atom pairs in an optical lat- 
tice even exemplifies stable states without any analogue 
in traditional condensed matter physics [TH - 

So far, a number of schemes have been proposed to de- 
rive an effective Hamiltonian to describe the spin-related 
dynamics for the multi-component system in the Mott- 
insulator phase 0, [1]. Most of the schemes ignore the 
existence of the upper bands and take single-band ap- 
proximation, which is reasonable when the on-site inter- 
action is much smaller than the energy gap between the 
first band and the second one. The situation may change 
dramatically if the scattering strength of the atoms is 
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greatly enhanced by the Feshbach resonance so that the 
on-site interaction exceeds the band gap. Recently, Kohl 
et. al. have studied the fermionic mixture of two hyper- 
fine states of in a three-dimensional optical lattice 
and accessed the strongly interacting regime via a Fes- 
hbach resonance, in which coupling between the lowest 
energy bands was dynamically generated Theoreti- 
cally, Diener and T.-L. Ho analyzed that a band insulator 
may evolve into the state with more bands occupied near 
the Feshbach resonance [I^. Very recently, A. F. Ho 
studied the phase transition from band insulator to Mott 
insulator for a fermionic system in optical lattices at a fill- 
ing of two fermions per site under the two-band approxi- 
mation . In that work, the Hund-like orbital coupling 
term is shown to play a special role in the strongly inter- 
acting regime and favors spin alignment between different 
orbits. 

It is thus physically nontrivial to go beyond the single- 
band approximation. Motivated by the recent progress 
on the research of the atomic gas in optical lattice near 
a Feshbach resonance, in this paper we study the equal- 
mixing two-component bosons in optical lattice with a 
filling of two bosons per site, focusing on the Mott- 
insulating regime and the spin-related phase transition 
due to the Feshbach resonance. As in the fermionic case, 
on each site there are many orbits and higher orbits may 
be occupied when the system is near the Feshbach reso- 
nance. Without loss of generality and for the purpose of 
simplicity, we take into account only two bands in the fol- 
lowing text, which can be fulfilled by enforcing the on-site 
interacting energy smaller than the energy level spacing 
between the third and the first orbital. We will show that 
in the strongly interacting regime the induced intcr-band 
coupling prefers the two atoms in different orbits on a site 
to form an on-site singlet, which is quite different from 
the Hund-like orbital coupling in the fermionic systems 
For simplicity, we consider only a one-dimensional 
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(ID) system which can be achieved by tuning the laser 
amphtudes Vn^^ ^ Vq^, Vq^ to produce a set of uncoupled 
ID tubes [3, llBl- In each tube, the system is effectively 
described by a ID optical lattice because the transverse 
motion is completely frozen. 

We start with the microscopic Hamiltonian of the two- 
component bosonic system in a ID optical lattice 

H = y2 r {x) d^ij^ (x)) 
+y2 (^0=^ ^^^^ ~ ^°') (^^^ 

+ Y.-aI '^^^i (^) "ti' (^) (^) (^) ' (1) 

<T,cr' "'0 

where the "spin" indices a [ indicate the two species 
of atoms or, equivalently, atoms with two internal states 
and /icr is the chemical potential. For the equal-mixing 
bosons, we have fii ^ fii ~ fj, corresponding to iV| = , 
where N^- is the total atom number of each specie. The 
optical lattice potential has the form of Vox sin^ kx with 
wave vectors k = 2tt/ \ and A the wavelength of the 
laser light. The parameter c > describes the repulsive 
interaction of the atoms and the interaction strengthes 
of intra-species and of inter-species are taken to be the 
same. Since we are interested in the regime where the 
interaction energy is tuned so that at most two Bloch 
bands are populated, it is sufficient to expand the oper- 
ator (x) in the lowest two Wannier functions 



ipa- (x) 



i,Q=1.2 



(2) 



where the operator Ci^a annihilates an atom with spin 
a in the band a at lattice site i. In a deep lattice the 
Wannier functions tOia (x) can be approximated by the 
local harmonic oscillator orbits in the ground state and 
the first excited state 
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, . (-1)' ^f2{x~x,) 
\x) = - — — ^ exp ■ 



(vra^) 



{x - XiY 



(3) 



(4) 



where oq = y^h/mujT is the ground state size of the local 
harmonic oscillator. Here ujt = V^yoxEa/h and Ej^ = 
ti^k'^/2m is the recoil energy. 

The second quantized Hamiltonian thus consists of 
three parts 



H ^ Ht + Hintra + Hi- 



(5) 



The hopping term Ht describes tunnelling of atoms from 
one site to another, which is typically assumed to occur 
between the nearest neighboring sites 



Ht = - 



E 



H.c. 



(6) 



and the hopping energy is 
ft2 rL 



tap 



2m 



dxdxUJi+ia (x) dxUJip (x) . (7) 



Hintra IS the contact-typc interaction Hamiltonian in the 
same energy band 



i.a.a i,a 

where the chemical potentials for each band 

r-i / fc2 



(8) 



are distinguished by a difference A — fii ^ fi2- This dif- 
ference is roughly the band gap between the two bands 
for deep lattice. On the other hand, the on-site interac- 
tion Hamiltonian between the two bands is denoted as 

Hinter 



+ 1 J2 Uo.p(n,ac.n,,f} + A'LA[.p)- (10) 

where A^^ = CiipCi^p, A^^^ = CiapCiap and 5*+ = 

'^l^a^iia {s^^ = is a pseudo-spin operator. The 

repulsive interaction (positive scattering length) between 
two atoms sharing a lattice site in the same band or be- 
tween the two bands gives rise to an interaction energy 



dxuiL ix) 



or 



(11) 



(12) 



which is just the additional energy that one needs to 
put two atoms on one site, in the same band or in dif- 
ferent bands. The term of S^S~jj describes the orbital 
coupling between the upper and lower bands/orbits. A 
striking feature here is that we have got an interaction 
term with opposite sign compared to the Fermionic case 
[l^ , for which the Hund-like orbital coupling term favors 
the "spin" of the two fermions at each site residing in 
different bands aligning paralleled. The orbital coupling 
thus determines the ground state in a different way for 
the bosonic case. The spins tend to align anti-paralleled 
in different bands when the interaction exceeds the en- 
ergy gap far away as illustrated later in Figure 2. The 
terms of Al^Aip and A^j^^A^^^ describe the interaction 
of atomic pair in different bands. 
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Substituting the approximate Wannier functions Eqs. 
© and (gl) into ^ and (HH) we easily obtain the 
parameters Un = U, U22 = 0.75U, U12 = 0.5U where 
U = c/4-\/27rao. Unlike the long-range Coulomb inter- 
actions for electrons in solids, here the orbital coupling 
term is of the same order of magnitude as the on-site 
repulsion term. Owing to the approximation of local 
harmonic oscillator orbits on Wannier functions, the in- 
tegral of the hopping matrix element between different 
bands is nonzero and in fact they satisfy the relations 
|iii| < 1*12! < 1^221- In optical lattice, both the hopping 
term tap and the on-site interaction U depend on the am- 
plitude Vq of the laser field. In this work, we will focus on 
the Mott phase with a larger ratio of C// tap and study 
the ground-state phase transition due to the change of 
the on-site interaction. In principle, via the Feshbach 
resonance, one could tune the strength of interaction so 
that U < A or A < U < 2 A. In the former case two 
bosons occupy the lowest Bloch band while in the latter 
case one of the atoms in the lowest band would be forced 
into the higher excited band. 

In the strong coupling limit with tap <C Uap,A, it is 
instructive to first consider the on-site local Hamiltonian 
with tap — and then treat the hopping term ^ as 
perturbation. It is easy to diagonalize the local on-site 
Hamiltonian (both intra- and inter-band parts ([8|) and 
(fTOp with two bosons per site. The local spectrum are 
given by: 



£1 = 


-2/ii 


+ 


7U 


£2 = 




+ A, 




£3 = 


-2/ii 


+ 


U 
2"' 


£4 = 




+ 


-u, 


£5 = 


-2/ii 


+ A4 


7U 





Among the ten eigenstates, those corresponding to 
eigenenergy ei are three-fold degenerate and are given 
by 

1+) = -e|TT,o)-i-/|o,TT), 
|o> = -e|Ti,o) + /|o,iT>. 

h> - -eUi,0)+/|0,ii>; 

The state corresponding to £2 is a local singlet forming 
by the atoms in the upper and lower orbits 

N) = ^(IT,i)-U,T)); 

The states corresponding to £3 are 

\t+) = IT,T>, 
\t-) = 
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FIG. 1: (Color online) The on-site energies Si versus U/A 
where we have made a total energy shift of — 2/^i. 



and state corresponding to £4 is 



1 



|to) = -^(|T,i) + U,T)); 

Finally the states corresponding to £5 are again three-fold 
degenerate 

1+') = /|TT,o> + e|o,TT), 

|0') - /|Ti,0>+e|0,iT>, 
I-') = /Ui,0>+e|0,ii>. 

Here we have used the notation for the representation 
of eigenstates that the left of comma in the right bar is 
for band 1 and the right of comma is for band 2. For 

example, IttiO) = (''Ti) 1*^^ represents two atoms 

with spin of 1 in the lower orbit and |t,i) — cj]^cj2 |0) 
represents an atom with spin of | in the lower orbit and 
an atom with spin of | in the upper orbit. The coefficients 
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fulfill (? + p — \ and (? and describe the probability 
two atoms simultaneously stay at the lowest band and 
the upper band respectively. We have < < 0.0659 



for < (7 < A. When U / /S ^ 0, p ^ and thus 
the system goes back to the single band model. To give 
concrete examples, we note — 0.0006 for U/A = 0.1 
and p — 0.0169 for U/ A — 0.5. Hence in the weakly in- 
teracting regime the two atoms mainly stay in the lowest 
band. 

In Figure 1, we display the five eigenenergies as a func- 
tion of U/A. To get the phase diagram for the ground 
state, it is sufficient to identify the lowest two levels while 
£3, £4 and £5 always correspond to the higher bands. The 
competition of the lowest two levels gives rise to com- 
pletely different ground state structure of the system and 
the transition point Uc/ A ~ 1.19 is approximately de- 
termined by the energy level crossing of £1 and £2 . For 
U > Uc, the local ground state on each site is a singlet 
state \s) with the spins of the two bosons aligned anti- 
paralleled. For U < Uc, the local ground state is one of 
the spin triplet \+), |0) and |— ). It is worthwhile to indi- 
cate that although the total spin fulfils (5*°*"') = as a 
result of — ^i, at each site the two species of bosons 
are not necessarily equal-mixing. At the first sight, this 
seems to imply that, in the limit of tap = and U < Uc, 
the ground state of the whole system is highly degener- 
ate and the spins of atoms at each site align arbitrarily 
because the local ground state on an isolated site can be 
either of the three states as long as the total spin of the 
system is zero. Actually this is not true when the hopping 
processes between the neighboring sites are considered. 

Now we switch on the hopping term between the 
nearest-neighbor sites. For the system with a filling fac- 
tor two, the state with two atoms at each site has low- 
est on-site energy. The process of an atom hopping to 
its neighboring site would change the on-site population, 
however such a hopping process is greatly suppressed be- 
cause placing three atoms at a site extremely costs en- 
ergy. Nevertheless, the virtual process of hopping to an 
intermediate state and then hopping back gives a second 
order correction to the ground state energy and lowers 
the ground state energy. The virtual hopping process 
does not change the total on-site populations but can ex- 
change two different atoms on neighboring sites. These 
virtual exchanging processes can be described by an ef- 
fective Hamiltonian acting on the ground states which is 
obtained in a second-order perturbation theory as 



Sr^ {^J-u+i\Ht \m) (m| Ht , w , 

2^ 7^ — 7^ {vi,i+i\ 

i.m 



Eq — Er, 



(13) 



, = \g)^®\g)^+i} are ground states 

1 = 2 and with ground state energy Eq. 



where {\^Ji^,i+l 
with Ni = N, 

As U < Uc, the three- fold degenerate ground states at 
an isolated site form a triplet, i.e., \g)^ = , |0)j , or 
|— )j, and therefore is nine-fold degenerate with 

Eq = 2e\. On the opposite regime oiU > Uc, = 
k)i '^it'^ ^0 = 2£2. The intermediate states are 

the product of states on the two neighboring sites with 
three and one bosons respectively and with excitation 
energies Em. 




-> JJ 



FIG. 2: Schematic picture for competing ground states in op- 
tical lattice. A phase transition from spin exchange to bosonic 
singlet occurs at U — Uc 



The second order perturbation calculation of the hop- 
ping terms enables us to identify one type of spin related 
quantum phase transition induced by the Feshbach res- 
onance. On the one side of the transition point, that is 
in the weakly interacting regime ([/ < Uc), the effective 
Hamiltonian can be further simplified and represented as 
an effective isotropic Heisenberg model in terms of spin-1 
operators. After straightforward but tedious algebra we 
get the effective Hamiltonian 



(id) 



(14) 



where S°' is a spin-1 operator in a (a = x,y,z) orienta- 
tion and the spin exchange coefficient is 



2 f2 



U 



2e^f 



A + U 3A- 



\U 



(15) 



where terms is neglected due to its smallness. Eq. 
(fT4|) is nothing but the Hamiltonian of an isotropic 5 = 1 
ferromagnetic quantum Heisenberg spin system. In the 
limiting case A ^ U , that is, when the upper band lies 
much higher than the lower band, we find 



A = 2 



u 



(16) 



In this limit, the probability of two atoms occupying the 
lowest band approaches unity. We then recover the re- 
sult in the single band approximation ■ The isotropic 
ferromagnetic model ([Mil has (25*°*°' -|- l)-fold degener- 
acy with 5*°*"' = A^/2 (A^/2 = 7V| = iV;). The ground 
state corresponds to the state with 5*°*"' = 0, in which 
case no spatial broken symmetry occurs. In the bosonic 
langauge, this means that the system phase does not sep- 
arate in the ground state for the SU (2)-symmetric model. 

On the other side of the transition point, when U > Uc, 
the ground state at the isolated site is a singlet. In this 
case, the virtual hopping process does not induce redis- 
tribution of on-site spins and the global ground state is 
the product of on-site singlets. We straightforwardly ob- 
tain the correction to ground state energy per site of the 
optical lattice by calculating the virtual hopping process 
to second order perturbation 



e = 62 + Se 



(17) 
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with 

2{2U ^ lU A + lU^ 2U-AJ- 

(18) 

Obviously this correction is negative and the hopping 
process always lowers the ground state energy. 

Figure 2 depicts the phase diagram of two-boson in 
two-band optical lattice model. For the interaction 
U < Uc, the atoms on a site form a triplet and the vir- 
tual hopping process produces ferromagnetic exchange 
between spins on neighboring sites, while in the strong 
coupling limit U > Uc the atoms in different bands prefer 
to align their spin anti-paralleled and form an on-site sin- 
glet. A phase transition from spin exchange to bosonic 
singlet occurs therefore dX U — Uc- We recall that in 
the fermionic case the phase diagram exhibits drastically 
different structure. Fermions with a filling factor two in 
two-band optical lattice are shown to exhibit opposite be- 
havior and there exists a phase transition from the band 
insulator to a Mott insulator with interesting dynamics 
of a spin-1 Heisenberg anti-ferromagnet p^ . 

We notice here a big difference between ID fermions 
and bosons. According to Haldane's conjecture, the 
ground state of 5'?7(2)-symmetric anti-ferromagnetic 
spin-1 Hamiltonian is gapped. Hence, small deviations of 
the Hamiltonian parameters reducing the SU (2) symme- 
try ioU{l) will not lead to qualitatively different results. 
On the other hand, the ground state for the ferromagnetic 
Hamiltonian is ordered and it is crucially important to 
consider a generic non-S'C/ (2)-symmetric model in order 
to understand whether the ground state is an easy-axis 
or easy-plane ferromagnet. 

To do this, we let the tunnelling matrix elements t^^ 
(ct —] , I) depend not only on band indices a, (3 but also 
on the component index a. Furthermore we distinguish 
the intra-species interaction U — U^^ = U^ and inter- 
species interaction U' = U^i to break the SU{2) sym- 
metry. When the system is in the strongly interacting 
regime, deviation of the SU{2) symmetry does not lead to 
qualitative change of the ground state properties because 
the ground state is composed of on-site singlets. How- 
ever, in the weakly interacting regime, when the SU (2) 
symmetry is broken, the effective Hamiltonian can be of 
the easy-axis type or of the easy-plane type with differ- 
ent kinds of ground states. We note that, for the general 
case with U ^ U\ the effective Hamiltonian can not be 
represented in the form of a simple spin exchange model. 
However, if |[/' — [/| ^ C/, C/', we can attribute the differ- 
ence of the on-site interacting energies to the zeroth-order 
Hamiltonian [H] and get an effective Hamiltonian of XXZ 
model 

+bY,S,, + dY,{S,,)\ (19) 

i i 



where 




In the limiting case A ^ [/, it is easy to show 
that A' = 2t{^t\^/U, SX'^ = [tl^ - tiilVC/, B = 
-3[(iii|)2 - itiiif]/U and D ^ U - U' and we recover 
the result for the single band approximation 5, 8]. It is 
obvious that we have always a positive small anisotropy 
parameter S\'^ for t]^^ ^ t''^^, which implies the effective 
XXZ model describes an easy-axis ferromagnet. Under 
the condition of 6X'^ ^ D, the ground state of the spin 
system is in a phase with spin domains. In the bosonic 
language, it corresponds to the situation with phase sep- 
aration of the two components. This implies that differ- 
entiating the tunneling terms for different components 
would induce phase segregation. When = we 
have SX'^ — and B = which reduces the model to 
(|14p except an additional term D (which vanishes natu- 
rally for S't/(2)-symmetric model because U — U'). For 
a large positive D, however, an easy-plane ground state 
can be realized. In terms of the nomenclature in Ref. 

, the easy-plane ferromagnet means the counter-flow 
superfluid. Straightforwardly, a positive D reduces the 
Sz component of the spin on each site. At large enough 
-D > 0, all spins will be essentially confined to the state 
with {Siz) — 0, which implies that large enough intra- 
component interaction (C7 ^ U') leads to two atoms 
belonging to distinct species occupying each site. On 
the other hand, for small enough D < Q {U <^ U'), the 
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ground state would stay in the state with {Siz) = ±1 and 
the term of D enhances the phase separation of different 
components. 

Before ending the discussion, we would like to remark 
the extension of the present work to the case with higher 
dimensions. Unlike the single band model which can be 
directly extended to the high-dimensional case, the effec- 
tive Hamiltonians (|14p and are no longer applicable 
to the high-dimensional optical lattice models when the 
higher orbits are populated. For higher dimensions, the 
first excited state in a local site is degenerate and has 
spatial anisotropy. Correspondingly, the hopping ma- 
trix element acquires spatial anisotropy and new phys- 
ical phenomena may arise due to the orbital degeneracy 

In summary, we have studied the quantum phase tran- 
sition induced by effective orbital coupling in optical lat- 
tices for an equal-mixing two-component boson system at 



a filling factor two per site. In the regime with weak on- 
site interaction, the two atoms stay in the lowest band 
and can be described by an effective spin- 1 ferromag- 
netic exchange model. In the regime with strong on-site 
interaction, the two atoms prefer to occupy different or- 
bits on a site and form an on-site singlet due to the ef- 
fective orbital coupling. We also considered the generic 
non-S'[/(2)-symmetry model. In the weakly interacting 
regime, the ground state may be described by an easy- 
axis ferromagnet corresponding to the case of phase sep- 
aration or an easy-plane ferromagnet corresponding to 
the state of counter-flow superfluid. 
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